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Abstract 

We study the quantum Hamilton- Jacobi (QHJ) equation of the recently obtained exactly 
solvable models, related to the newly discovered exceptional polynomials and show that 
the QHJ formalism reproduces the exact eigenvalues and the eigenfunctions. The fact 
that the eigenfunctions have zeros and poles in complex locations leads to an unconven- 
tional singularity structure of the quantum momentum function p(x), the logarithmic 
derivative of the wave function, which forms the crux of the QHJ approach to quantiza- 
tion. A comparison of the singularity structure for these systems with the known exactly 
solvable and quasi-exactly solvable models reveals interesting differences. We find that 
the singularity structure of the momentum function for these new potentials lies between 
the above two distinct models, sharing similarities with both of them. This prompted us 
to examine the exactness of the supersymmetric WKB (SWKB) quantization condition. 
The interesting singularity structure of p(x) and of the superpotential for these models 
has important consequences for the SWKB rule and in our proof of its exactness for these 
quantal systems. 

1 Introduction 

New infinite sets of solvable quantum mechanical potentials have been recently con- 
structed by Odake and Sasaki PQ, [2], using the Darboux-Crum transformations [3], [I]. 
They deformed the standard radial oscillator and Darboux-Poschl- Teller potentials using 
the Laguerre and the Jacobi polynomial eigenfunctions of degree I. The solutions of the 
new potentials are in terms of the newly discovered Laguerre- and the Jacobi- Xi type 
exceptional orthogonal polynomials respectively, which have been recently introduced by 
Gomez-Ullate et.ai, [5], [6]. The characteristic feature of this new class of polynomials 
is that the series starts with a polynomial of degree £(£=1,2...) and still can form a 
complete set with respect to a positive-definite measure. This is unlike the classical or- 
thogonal polynomial series, which require a constant to form a complete set. In addition 
to this feature, the eigenfunctions of these new potentials have singularities at complex 
locations, associated with the zeros of the eigenpolynomials of the standard potentials 
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from which they are derived. These standard polynomials seem to play a crucial role 
in the construction of these new potentials. The fact that the new Laguerre- and the 
Jacobi- type polynomials occur as solutions to the Sturm-Liouville problem with ra- 
tional coefficients, associated to the standard Laguerre and Jacobi polynomials, is very 
intriguing. In this light the quantum Hamilton- Jacobi (QHJ) analysis [7J, [8], which 
uses the singularity structure of the logarithmic derivative of the wave function to obtain 
the eigenvalues and eigenfunctions for a given potential, can reveal interesting features 
pertaining to these models. 

In this paper, we analyze the set of potentials obtained by deforming the radial oscilla- 
tor using the quantum Hamilton- Jacobi (QHJ) formalism and investigate the singularity 
structure of the logarithmic derivative of the wave function, i[)(x): 



known as the quantum momentum function (QMF). The knowledge of the singularity 
structure of the QMF, coupled with the exact quantization condition, 



defined within this formalism, allows us to arrive at the required solutions [9]- [12]. Here, 
n gives the number of nodes of ip(x) and the contour C encloses these nodes located in the 
classical region, in between the two turning points x\ and 22 in the complex x-plane. In 
the present study, in addition to obtaining the energy eigenvalues and eigenfunctions, we 
concentrate on bringing out the new features of the singularity structure of the QMF of 
these potentials. This, when compared to that of the models studied earlier, namely the 
exactly solvable (ES) models, whose solutions are in terms of the classical polynomials, 
and the quasi-exactly solvable (QES) models revealed interesting differences. This led us 
to a careful analysis of the exactness or otherwise of the supersymmetric WKB (SWKB) 
quantization condition [13], [2] for these new models, because it was shown in [15] that 
the singularity structure of the QMF provided a link to the exactness of the SWKB 
integral. Here, Bhalla et.al., showed that the SWKB condition is exact for all models, 
where the locations and the residues of the poles in the non-classical regions of the SWKB 
integrand and the QMF matched identically. Hence, it will be interesting to see if a similar 
analysis works for these new models, especially when their singularity structure differs 
from that of the conventional models. 

In the following section, we give a brief account of the QHJ formalism [7]- [12], fol- 
lowed by the analysis of the generalized deformed radial oscillator potential within this 
method. We then compare the singularity structures of various models, in order to iden- 
tify the similarities or differences, between the ES and QES models and the present ones. 
In section 3 we investigate the exactness of the SWKB condition for the deformed ra- 
dial oscillators and subsequently conclude in the last section, with directions for further 
research. 




x)dx = nh, 



(2) 
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2 QHJ formalism and the deformed radial oscillators 



In the QHJ formalism, x is treated as a complex variable and we use the Riccati equation 
for q (= ip(x)) (H = 2m = 1) 



q 2 + d x q + (E-V(x)) = 0, (3) 

<//(*) 

tp{x) 



to obtain the required solutions. Using (JTJ and substituting q = ^4^- in the above 



equation gives us the Schrddinger equation 

- ^-P^- + V(x)^(x) = EiP(x), (4) 
ax 1 

which establishes the connection between QHJ approach and the conventional ones. In 
this approach, the singularities of the momentum function q play the key role in obtaining 
the spectra, without the need to solve for the eigenfunctions. 

The singularities of q are of two types namely, the fixed and the moving singularities. 
The fixed singularities correspond to the singularities of the potential and their location 
is independent of the energy and can be located by inspection from flS]). The moving 
singularities are poles, which correspond to the zeros of the wave function ip(x). It is well 
known that the location of the nodes of the wave function changes with energy. Therefore, 
the location of these poles depends on the initial condition and hence they are named 
moving poles. The number and the location of these poles cannot be inferred from the 
differential equation by inspection. These poles turn out to be finite in number for all the 
ES models studied [10], [T2]. This is equivalent to the fact that the point at infinity, at 
best, is an isolated singularity. The location of all these singularities and their residues 
are used to evaluate the integral in fl2]), which gives the expression for the eigenvalues for 
the ES models [10] and the quasi-exact solvability condition for the QES models [TT] . In 
addition, the knowledge of the singularity structure of q allows us to obtain the form of 
QMF and hence the eigenfunctions using ([T]). For more details, we refer the reader to 
our earlier papers [H]-[I2] and we proceed to analyze the potentials under study. 

Deformed radial oscillator potentials 

Here, we present a brief account of the set containing infinite number of new shape invari- 
ant potentials (SIPs) Vg(x), indexed by £ (£ = 1, 2, ....) described in [TJ, [2J. The potentials 
Vg(x) constitute the hierarchy of supersymmetric potentials and are translationally shape 
invariant [13], [16]. Moreover, within each set, putting £ = gives the respective standard 
potential and £ = 1 corresponds to the potentials constructed by Quesne [T7]. In [18], 
all the SIPs including the newly deformed potentials have been obtained from a general 
analysis of shape invariance requirement. Here, we analyze the set of potentials obtained 
by deforming the standard radial oscillator. For every potential in this set, indexed by £, 
a prepotential 

u e {x; g) = u {x; g + £) + In — [ X ■ 9 + ^ , g > (5) 



is defined, where 



coo(x; g) = —-x 2 + g\ogx, < x < oo, (6) 

Ux 2 ;g) = L^ l \-x 2 ). (7) 
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Here L e 2 (x) is the associated Laguerre polynomial. The corresponding superpoten- 
tial is given by 

ur f \ a / \ 9 + Z d x &(x 2 ;g + l) d x &{x 2 ;g) 

{x) = - & * W ^ 9)=x ~ — - U^,9 + i) + ~U^gJ (8) 

and the corresponding potential, Ve(x), with zero ground state energy is given by 

Vt(x) = (d x w e (x; g)) 2 + d xx w e (x; g). (9) 
Substituting ue(x; g) in the above equation and simplifying we obtain, 

V e (x) = U 2 (x) + d x U (x) + d: -' {J ' J ' : " : U 



+2 (-x + i±i - d fff ]9 A d fff ;9 + ]\ (10) 



where 

Tr , x , g + ? d x &(x 2 ;g) 

From the differential equation and the recurrence relations of the associated Laguerre 
polynomials [12], [2], we obtain 

W; | , + 1) = _ 2 ( + t _+t\ aMf;9 + i) +4t (12) 

&{x 2 ;g + l) \ x J &(x 2 ;g + l) 

and 

( dM(x 2 ;g) \ ( dM(x 2 ;g + l) \ = d x U^g + i) ( _ ^(x 2 ;^) \ , . 

U^ 2 ;<?) A ^ 2 ;5 + i) J x a* 2 ;g) V ^ 2 ;5 + i)J [) 

respectively The use of these equations in (fTOl . reduces the potential to 



V e (x) = U 2 (x) + d x U (x) - Ax dM }fl 9 | 1} + 4*. (I I) 

Substituting (IT41) in ([3]) gives us the QHJ equation 

q 2 + d x q t + E = U 2 (x) + d x U (x) - Ax U ff \ 9 ± 1} + At (15) 

Here, qi corresponds to the momentum function associated to Vi(x). From the above 
equation we can see that qe has a simple pole at the origin along with 2£ fixed poles 
corresponding to the zeros of ^(x 2 ;g). An important feature of all known ES models 
studied has been that the point at infinity is an isolated singularity [TD] . The consequence 
of this is that the QMF has finite number of moving poles. For the new ES models under 
study we make the same assumption, i.e., momentum function has a finite number, N, 
moving poles. From (TT5]) . we can see that qi behaves like x at infinity. This coupled 
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with the knowledge of its singularities and residues allows us to write qe in the form of a 
meromorphic function 



Qe 



-x + 



9 + 



x 



21 

E 



X 



N 



J= i * - bj 



(16) 



The first term on the right hand side of the above equation gives the large x behaviour 
of qi and the second term is the singular part arising from the pole at the origin. The 
numerator is the residue at the origin and is calculated by substituting the Laurent 
expansion of qe around the same : 



qe = h a + a x x + 

x 



;i7) 



in (TH)]) . Equating the coefficients of the same powers of x to zero, we obtain two values 
for b\, namely ±(g + £). The correct value is chosen using the boundary condition, in the 
limit E — > 0, qe — )> — We(x). This gives the residue to be (g + 1). The third term in ([16]) 
gives the singular part associated with 2£ fixed poles, with residue —1, corresponding to 
the zeros of the £e(x 2 ,g). It is easily seen that this is equivalent to g f ■ The fourth 
term corresponds to N simple moving poles with residue 1. Some of these are located 
on the real line and others off the real line. The ones which lie on the positive real line 
correspond to the nodes of the wave function. The contribution of the moving poles of 
qe is similarly written as %^~t~q > with Pn(x) being a polynomial of degree N. The last 
term, <p(x), describes the analytic part of qe. From Liouville's theorem, <p{x) is a constant, 
which turns out to be zero from (Tl5l) . Substituting (fl6|) in (fl5j) and after a few algebraic 
manipulations we get 

d xx P N (x) 2 (_ x + g + £ d x &(x 2 ; g)\ d x P N (x) 



Pn(x) 



+4x 



x ie{x 2 ;g) J P N (x) 
dxie{x 2 ]g + l) 



&(x 2 ;g) 



M+E = 0. 



(18) 



For large x the leading term in the above equation is a constant, therefore equating it to 
zero, we obtain 

E = 2N — At (19) 



Substituting this in ([18]) . we obtain 



d xx P N (x) + 2 



-x + 



9 + £ _ dx£,i{x 2 -g) 
x ie(x 2 ;g) 



&(x 2 ;g) 

The above equation can be written as 



[x dM* 2 - : 9 + i) +2N _ 8i ) Pn {x 



d x P N (x) + 
0. 



(20) 



1 d 



d 



2 j \ w e(x)—P N (x) 
wf dx \ dx 



Ux 2 ;g) J 



(21) 
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where, 



WAX) 



cxp 



exp 




dx€e(x 2 ;gY 



(22) 



The differential equation (12 ip implies that wf is the weight function with respect to which 
the polynomials, Pn{x), are orthogonal. The function wf (x), coincides with the weight 
function, associated with the exceptional Laguerre polynomials [1]. Therefore, apart 
from an overall constant, Pn{x) are the exceptional Laguerre polynomials p£ tn (x 2 ;g) 
with degree 2n + 21 in x pp. This implies N = 2(n + £) , which when substituted in (Tl9|) 
gives the energy eigenvalue expression to be 

E = An. (23) 

With N = 2(n + £), it is easy to verify that (f20|) . coincides with the differential equation 
[2] for the exceptional Xi Laguerre polynomials 

d x &(x 2 ;gY 



d xx P N (x) +2[-x + 



g + 



An-A£ + Ax 



x &(x 2 ;g) 
dM{x 2 ,g + l) \ 
Ux 2 ,g) J 



d x P N (x) 



Pn(x) 



0. 



(24) 



In addition, substituting the meromorphic form of the momentum function from (I16p in 
(PQ), we obtain the unnormalized wave function as 



x 



x 



&(x 2 ;g) 



9), 



(25) 



which agrees with the known results pQ. Thus, we have obtained the energy eigenvalues 
and eigenfunctions for the new exactly solvable potentials, using the singularity structure 
of the momentum function. In addition to the fixed poles, qe has 2n moving poles on the 
real line and 21 moving poles off the real line. The point at infinity turns out to be an 
isolated singularity. 



Singularity structure 



Here, we compare the singularity structure of the momentum function associated 
with these models to those of the QES and the conventional ES models. We see that 
the differences among these three models are due to the moving pole structures. For all 
the models, the number of moving poles located on the real line is consistent with the 
oscillation theorem. For the conventional ES models there were no moving poles off the 
real line. For the QES models, moving poles off the real line were found. The number of 
these poles varied from one QES level to the other, keeping the total number of moving 
poles (located at both real and complex locations) fixed for all solvable states. Compared 
to this, the new ES models have moving poles off the real line, but their number is 
fixed, equal to 2£, for all the energy levels. Thus, we see that the singularity structure of 
the newly constructed potentials is neither completely like the ES models nor the QES 
models. 
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3 Exactness of SWKB quantization condition 



In an earlier paper [15], the information on the singularity structure of the QMF for several 
models, was used to account for the exactness of the SWKB quantization condition, 

- r J[E- W 2 (x)]dx = n, (26) 

TT J xi 

where x\ and x 2 [x\ < x 2 ) are the turning points, which are the real roots of Je — W 2 (x). 

For simplicity we denote the integrand \J E — W 2 {x) as Pswkb and put h = 1. The 
integrals such as fl26|) can be turned into a contour integral (see for example [20]), and 
for fl26|) we get 

l c = 7T f VswKEtdx = n, (27) 
2tt Jc 

where C is the contour enclosing the branch cut of Pswkb from x\ to x 2 . This integral 
can be computed in terms of the contributions of the singularities outside the contour 
C. The residues of Pswkb are double valued and the right values are picked using the 
boundary condition, in the limit E — Y 0, Pswkb — > P, since W(x) = — tp' (x) / ipo(x) . 
For the cases studied earlier in [15], the residues and the locations of the singularities 
outside the contours C of the two integrals, f[2T|) and ([2]) matched identically and hence 
application of complex integration techniques gave same answers. 

In the present case, we first list out the location of the poles of pswkb outside the 
contour C and corresponding residues in table 1 below. It can be easily seen that the 



location 


residue 


x = 


-i(g + £) 


21 singularities 


—i 


located at 




the zeros of 




&(x 2 ,<7 + l) 




2£ singularities 


% 


located at 




the zeros of 








x = oo 


f + i(g + e) 



Table 1: Pole structure of Pswkb and p(x). 



poles and their residues are identical to those of the QMF. However in this case, Pswkb 
has additional branch points as illustrated in figure 1, for the case of I — 1 and g — 1. 
In figure 1, the contour 9 encloses all the poles and branch cuts of Pswkb- The con- 
tour C\ encloses the origin and the contour C 2 encloses the branch cut between — X\ and 
— X2, which is equal to the integral in (12~T|) by symmetry. The contours j m , T m enclose 
the simple poles corresponding to the 2£ zeros of £,e(x 2 ,g) and ^(x 2 ,g + 1) respectively. 
Finally, the contours Qi enclose the branch cuts of Pswkb located off the real line. By 
writing the left hand side of (!2~T|) in terms of contour integrals enclosing the poles and 
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Figure 1: The figure displays the singularities and branch cuts of J E — W 2 (x), for g = 1 



and I = 1 and the contours enclosing them. All the contours here are oriented in the 
anti-clockwise direction. 

branch cuts of Pswkb located outside C gives 



We show that the contribution of the branch cuts of pswkb enclosed by contours Q{ 
should vanish, for the SWKB condition to be exact. For this purpose we will evaluate 
the different integrals, except the ones around the contours Qi in the left hand side of 
( I2"gj) . The integral around contour is equal to 2ni( residue at x equal to oo), calculated 
using the rule [2U] 




(28) 



8 



Residue of J E — W 2 (x) at z equal to oo = — Coefficient of t in the 



expansion of yE — W 2 (l/t) for small t. (29) 
BY using f^±i) = E?£i _1_ and ^M|M = E f =t _i_ one obtains 



" t * t ' 2 



E-W 2 m ^- t ^-^-( 9 + ^-t^ — + t^—j , (30) 



21 f2 21 f2 

l-E*-2 {g + - g — + £ — + *( t ), (31) 



.7 



where -ff(t) denotes all the other higher order terms. Expanding the square root, one 
obtains 

1 



E - W \l/t) =±-[l--( E + 2(g + £Y)t 2 + ...,) (32) 

which gives the residue at infinity to be ±i(^ + (g + £)) . The application of the boundary 
condition E — > 0, Pswkb P, allows us to choose the right value of the residue to be 
Kf + (9 + 1))- Therefore 



2tt 

'f + (* + *)) (33) 

The other method to calculate the integral around G is to use the mapping x = 1/t, which 
maps the point at infinity in the x plane to the origin in the t plane. The coefficient of 
1/t in the Laurent expansion of the integrand in powers of t, will be the residue and is 
same as the value calculated above. 

The residues at the other singular points, listed in table 1, can be obtained by using 
simple complex variable techniques. Now with all the residues known, application of the 
Cauchy residue theorem gives, 



i c = j + (g + t) 



({g + £) + I c + [~2£ + 2£] + W Psw KB dx) . (34) 



Simplifying the right hand side gives 



E 1 u f 

2I c = tt ~ ttY. ? PswKBdx. (35) 
2 2n r-[ 



If E coincides with the n eigenvalues given in f[2"3"j) . we obtain 

u 



Ic = n - f PswKBdx. (36) 
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Therefore, it is obvious that if the sum of the contribution of the branch cuts off the real 
line vanishs, we obtain 



Thus vanishing of the contribution of the branch cuts off the real line turns out to be a 
necessary and sufficient condition for the exactness of the SWKB rule. 

We point out here that the calculation of the integral appearing in in terms of p, 
also proceeds in a similar fashion. The only difference being that the QMF does not have 
branch points or moving branch point singularities in the complex plane. Therefore this 
integral receives contribution only from the fixed poles corresponding to the potential, 
the n moving poles located symmetrically between —x\ and — x<i and from the point at 
infinity. As the table suggests the contribution from these poles, for both integrals is 
identical. Therfore it is obvious that the SWKB condition is exact if the contribution 
from the branch cuts off the real line vanishes, when the energy is an eigenvalue. 

In [16], [21] it has been proved that the SWKB condition is exact for all the conven- 
tional ES, translationally SIPs. The proof makes use of the property of shape invariance, 
where the assumption that W 2 (x) is of 0(h°) while hd x W(x) is of 0(h) was the key 
point. However, unlike the conventional SIPs with translation, for the new SIPs with 
translation, W(x) is in fact a complicated function of h. Thus the derivation in [16], [21] 
about the exactness of SWKB for SIPs with translation is strictly not applicable for the 
newly discovered SIPs with translation. Our present study of the exactness of the SWKB 
rule suggests that the disappearance of the contribution of the branch cuts off the real 
line, when E is equal to the one of the eigenvalue is possibly linked to the property of 
shape invariance. 

4 Conclusions 

In this study, we have analyzed the new ES shape invariant potentials obtained by de- 
forming the radial oscillator using the QHJ formalism and obtained the expressions for 
the eigenvalues and eigenfunctions. We have investigated the singularity structure of the 
QMF and compared it with those of the conventional ES and QES models studied earlier. 
The QMF for these potentials has 2n real and 2£ complex moving poles, with E = 4n 
giving the energy eigenvalue and £ fixing the potential. We found that the singularity 
structure of these models is neither completely like the ES models nor the QES models. 
The exactness of the SWKB rule has been shown to be equivalent to vanishing of the 
contribution of the branch cuts off the real line. Lastly, we point out that the above 
results can be extended to the other families of infinite number of potentials obtained by 
deforming the Darboux-Poschl- Teller potentials [TJ. 
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